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Abstract 

This paper is about cubic sand grains moving around on nicely packed 
columns in one dimension (the physical sand pile is two dimensional, but 
the support of sand columns is one dimensional) . The Kadanoff Sand Pile 
Model is a discrete dynamical system describing the evolution of a finite 
number of stacked grains — as they would fall from an hourglass — to a 
stable configuration. Grains move according to the repeated application of 
a simple local rule until reaching a stable configuration from which no rule 
can be applied, namely a fixed point. 

The main interest of the model relies in the difficulty of understanding 
its behavior, despite the simplicity of the rule. We are interested in de- 
scribing the shape of fixed point configurations according to the number of 
initially stacked sand grains. In this paper, we prove the emergence of a 
wavy shape on fixed points, i.e., a regular pattern is (nearly) periodically 
repeated on fixed points. Interestingly, the regular pattern does not cover 
the entire fixed point, but eventually emerges from a seemingly highly dis- 
ordered segment. Fortunately the relative size of the part of fixed points 
non-covered by the pattern repetition is asymptotically null. 
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1 Introduction 



1.1 Framework 

Understanding and proving regularity properties on discrete dynamical systems 
(DDS) is readily challenging, and inferring the global behavior of a DDS de- 
fined with local rules is at the heart of our comprehension of natural phenomena 
|Wea481 [HBLJ09J. There exists a lot of simply stated conjectures, often issued 
from simulations, which remain open (for example the famous Langton's Ant 
[GPST951 IGMG02] ). Sand pile models are a class of DDS defined by local rules 
describing how grains move in discrete space and time. We start from a finite 
number of stacked grains — in analogy with an hourglass — , and try to predict 
the asymptotic shape of stable configurations. 

Bak, Tang and Wiesenfeld introduced sand pile models as systems presenting 
self-organized criticality (SOC), a property of dynamical systems having critical 
points as attractors |BT88| . Informally, they considered the repeated addition 
of sand grains on a discretized flat surface. Each addition possibly triggers an 
avalanche, consisting in grains falling from column to column according to simple 
local rules, and after a while a heap of sand has formed. SOC is related to the 
fact that a single grain addition on a stabilized sand pile has a hardly predictable 
consequence on the system, on which fractal structures may emerge |Cre96| . This 
model can be naturally extended to any dimension. 

A one-dimensional sand pile configuration can be represented as a sequence 
(^i)ieN °f non-negative integers, hi being the number of sand grains stacked on 
column i. The evolution starts from the initial configuration h where ho = N 
and hi = for i > 0, and in the classical sand pile model a grain falls from 
column i to column % + 1 if and only if the height difference hi — hi + \ > 1. One- 
dimensional sand pile models were well studied in recent years [GK93, DL98, 
IGLM+04I [FMP071 lPha08l IPPP111 IFPPHllj . 



In this paper, we study a generalization of classical sand pile models, the 
Kadanoff Sand Pile Model. A fixed parameter, denoted p, is defined as the 
number of grains falling at each rule application. The results we expound in 
this paper present an interesting feature: we asymptotically completely describe 
the form of stable configurations, though there is a part of asymptotically null 
relative size which remains mysterious. Furthermore, proven regularities are 
directly issued from this messy part. This point and its link to the half-discrete 
half-continuous nature of sandpile models will be discussed in the conclusion. 



We formally define the model in subsection 1.2 and present a broader class of 



DDS in which it belongs in subsection 1.3 Subsection 1.4 exposes an interesting 



and useful way of computing fixed points, and subsection 1.5 introduces the 
results we expound in section [2] At the light of those developments, section [3] 
discusses that sand pile models exhibit a behavior at the edge between discrete 
and continuous phenomena. 

In some previous works ([PRllaJ[PRllb|), we obtained a similar result for 
the smallest parameter p = 2 (the case p = 1 is the well known Sand Pile Model), 
but, for the general case, we have to introduce a completely different approach. 
The mains ideas are the following: we first relate different representations of a 
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Figure 1: KSPM(p) rule. When p grains leave column i, the height difference 
bi-i is increased by p, bi is decreased by p + 1 and bi +p is increased by 1. Other 
columns' height differences are not affected. 



sand pile configuration (subsection 2.1 ), which leads to the construction of a DDS 
on Z p+1 such that the orbit of a well chosen point (according to the number of 
grains N) describes the fixed point configuration we want to characterize. This 
system is quasi-linear in the sense that the image of a point is obtained by a 
linear contracting transformation followed by a rounding (in order to remain in 
Z p+1 ). We want to prove that this system converges rapidly, but the rounding 
makes the analysis of the system very difficult (except for p = 2). The key-point 
(subsection 2.2) is the reduction of this system to another quasi-linear system 
on Z p , for which we have a clear intuition (subsection 2.3), and which allows to 
conclude (subsection 2.4). 



1.2 Definition of the Kadanoff Sand Pile Model 

Kadanoff et al proposed a generalization of classical models in which a fixed pa- 
rameter p denotes the number of grains falling at each step [KNWZ89] . Starting 
from the initial configuration composed of N stacked grains on column 0, we 
iterate the following rule: if the difference of height between column i and i + 1 
is greater than p, then p grains can fall from column i, and one grain reaches 
each of the p columns i + l,i + 2, . . . ,i + p. The rule is applied once during each 
time step. In the following, column and index are synonyms; and for the sake of 
imagery, we always consider indices (column numbers) to be increasing on the 
right (see figure [T]). 

Formally, it is defined on the space of ultimately null decreasing integer se- 
quences where each integer represents a column of stacked sand grains. Let 
h = (hi) j G Dsl denote a configuration of the model, hi is the number of grains on 
column i. 

In order to consider only the relative height between columns, we represent 
configurations as sequences of height differences b = (&i)i 6 N ; where for all i ^ 
0, bi = hi — hi+i. This latter is the main representation we are using, within the 
space of ultimately null non-negative integer sequences. 

We give two definitions of the model, obviously isomorphic. Definition [T] is 
more natural, but Definition [2] is more convenient and is the one we will use 
throughout the paper. 
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Definition 1. The Kadanoff sand 
pile model with parameter p > 0, 
KSPM(p), is defined by two sets: 

• configurations. Ultimately null 
non-negative and decreasing inte- 
ger sequences. 

• transition rules. We have a tran- 
sition from a configuration h to a 
configuration h! on column i, and 
we note h — ► h! when 

- h[ = hi-p 

- h' l+k = h i+k + 1 for < k «S 
P 

- h'j = hj for j f + 
l,...,i +p}. 



Definition 2. The Kadanoff sand 
pile model with parameter p > 0, 
KSPM(p), is defined by two sets: 

• configurations. Ultimately null 
non-negative integer sequences. 

• transition rules. We have a tran- 
sition from a configuration b to a 
configuration b' on column i, and 
we note b —> b' when 

- b'^ = h-!+p (for i # 0) 

- b^h-ip-1) 

~ t>i+p = h+p + 1 

- &'■ = bj for j $ {i-l,i,i+p}. 



We also say that i is fired. Remark Again, remark that according to the 
that according to the definition of the definition of the transition rules, i may 
transition rules, i may be fired if and be fired if and only if b . > p ^ otherwise 
only if hi - h i+1 > p, otherwise h[ is y. is negative, 
negative or the sequence h' is not de- 
creasing. 

We note b — > b' when there exists an integer i such that b ^* b'. The transitive 
closure of — > is denoted by — >, and we say that b' is reachable from b when b ^> V . 
A basic property of the KSPM model is the diamond property : if there exists 

i and j such that b b' and b ^> 6", then there exists a configuration b'" such 

that b' X b'" and 6" ^> 

We say that a configuration b is stable, or a /ixec? point, if no transition is 
possible from 6. As a consequence of the diamond property, one can easily check 
that, for each configuration b, there exists a unique stable configuration, denoted 
by vr(6), such that b vr(6). Moreover, for any configuration b' such that b ^> 6', 
we have 7r(6') = 7r(6) (see |GMP02] for details). We denote W the infinite se- 
quence of zeros which is necessary to explicitly write the value of a configuration. 
An example of evolution is pictured on figure [2] For convenience, we denote N 
the initial configuration (iV, W ), such that 7r(iV) is the height difference repre- 
sentation of the fixed point associated to the initial configuration composed of 
N stacked grains (see figure [3] for an illustration). This paper is devoted to the 
study of ir(N) according to N. 

A useful representation of a configuration reachable from (N, W ) is its shot 
vector (oi)ielMi where a, is the number of time that the rule has been applied 
on column i from the initial configuration (see figure [2] for an example) . This 
representation will play a major role in the following. 
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Figure 2: A possible evolution in KSPM(2) from the initial configuration for 
N = 24 to vr(24). tt(24) = (2, 1, 2, 1, 2, U ) and its shot vector is (8, 1, 2, (F). 



(24,0") 
I 

(21,0,1,0") 
I 

(18,0,2,0") 
I 

(15,0,3,0") 
I 

(12,0,4, 0") 
I 

(9,0,5,0" 
I 

(6,0,6,0" 
I 

(3,0,7,0") 
I 

(0,0,8,0") 
I 

(0,2,5,0,1,0") 



(15,2,0,0,1,0") 
I 

(12,2,1,0,1,0") 
I 

(9,2,2,0,1,0") 
I 

(6,2,3,0,1,0") 
1 

(3,2,4,0,1,0") 
I 

(3,4,1,0,2,0") 



(6,4,0,0,2,0") 
I 

,1,0,1,2,0") 



(0,4,2,0,2,0") (5,1,1,1,2,0") 



(2,1,2,1,2,0") 



Figure 3: The set of reachable configurations for p = 2 and N = 24. The initial 
configuration is on the top, and the unique fixed on the bottom. 
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Figure 4: The initial configuration of KSPM(2) is presented as a CFG where 
each vertex corresponds to a column (except the sink, vertices from left to right 
corresponds to columns 0, 1, 2, 3, ... ) with a load equal to the difference of height 
between column i and For example the vertex with load N is the difference 

of height between column (N grains) and column 1 (0 grain). 

1.3 Chip firing game 

Sand pile models are specializations of Chip Firing Games (CFG). A CFG is 
played on a directed graph in which each vertex v has a load l(v) and a threshold 
t(v) = deg + (v^ and the iteration rule is: if l(v) ^ t(v) then v gives one unit 
to each of its neighbors (we say v is fired). As a consequence, we inherit all 
properties of CFGs. 

Kadanoff sand pile is referred to as a linear chip firing game in [GMP02J. 
The authors show that the set of reachable configurations endowed with the 
order induced by the successor relation has a lattice structure, in particular it 
has a unique fixed point. Since the model is non-deterministic, they also prove 
strong convergence i.e. the number of iterations to reach the fixed point is the 
same whatever the evolution strategy is. The morphism from KSPM(2) to CFG 
is depicted on figure |4j 

When reasoning and writing formal developments about KSPM, it is conve- 
nient to think about its CFG representation where local rules let units of height 
differences move between columns. 

1.4 Hourglass, inductive computation and avalanches 

In order to compute tt(N), the basic procedure is to start from the initial config- 
uration (N, 0^) and perform all the possible transitions. However, it also possible 
to start from the configuration (0 W ), add one grain on column and perform all 
the possible transitions, leading to vr(l), then add another grain on column 
and perform all the possible transitions, leading to 7r(2), etc... And repeat this 
process until reaching ir(N). 

Formally, let 6 be a configuration, denotes the configuration obtained by 
adding one grain on column 0. In other words, if b = (bo, b\, . . .) then b^° = 
(bo + 1, bx, . . . ). The correctness of the process described above relies on the fact 
that 

(k,0 u ) A n(k- l) i0 

Indeed, there exists a sequence of firings, named a strategy, (si) l jZ} l such that 
(k — 1, W ) — * . . . ^* n(k — 1). It is obvious that using the same strategy we have 

deg + (v) denotes the out-degree of v. 
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(jfe-l,O w ) — nr(fc-l) 

(k, 0") > 7r(fe - l)+° > n{k) 



Figure 5: inductive computation of ir(k) from 7r(k — 1). The bold arrow represents 
an avalanche. 



nxr 



Figure 6: An example of avalanche: starting from vr(24), we add one grain on 
column (darkened on the leftmost configuration) and apply the iteration rule 
until reaching 7r(25). Arrows are labelled by the index of the fired column (the 
leftmost unstable column is fired at each step). 



(k, W ) = (k — 1, 0")''' — * . . . ^* n(k — 1)^° since we only drag one more grain on 
column along the evolution, which does not prevent any firing (see figure [5]). 
Thus, with the uniqueness of the fixed point reachable from (k, W ), we have the 
recurrence formula: 

7r(vr(A; - = ir(k) 

with the initial condition vr(0) = W , enabling an inductive computation of 7r(fe). 

The strategy from ir(k — 1)^° to Tr(k) is called an avalanche. Note that from 
the non-determinacy of the model, this strategy is not unique. To overcome 
this issue, it is natural to distinguish a particular one which we think is the 
simplest: the k th avalanche s k is the leftmost strategy from ir{k — 1)^° to 7r(fe), 
where leftmost is the minimal strategy according to the lexicographic order. This 
means that at each step, the leftmost possible firing is performed. A preliminary 
result of [PRllaj is that any column is fired at most once during an avalanche, 
which allows to write without ambiguity for an index i: i e s k or i ^ s k . 

A recent study ([GM10J) showed that in the two dimensional case the avalanche 
problem (given a configuration and two indices i and j, does adding one grain on 
column i has an influence on index j?) on KSPM is P-complete, which points 
out an inherently sequential behavior. 



1.5 Objectives of the paper 

Let us now introduce notations and state the results. Regarding regular expres- 
sions, let e denotes the empty word, and • (respectively +) denote the concate- 
nation (respectively or) operator. -\- is to the or what E is to the sur/j^] and 

k 

2 +0 ! =e + + 00+ - +0...0. 
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(a) We call wave the pattern 
p ■ . . . • 2 • 1 in a sand pile 
configuration. 




(b) Theorem [|] Illustration 
of the regular expression 
p+i \ 

+ •/"••••-'• I 

i=0 J 

as a sand pile configuration. 



0(ViV) 

(c) Theorem [Z] From an 
index in O(logiV) the fixed 
point 7r(JV) consists of waves 
all consecutive to each other, 
except at at most one place 
where two waves may be sep- 
arated by a 0. 



Figure 7: A wave (figure 7a), and graphical representations of Theorem [3] (figure 
7b ) and Theorem E] (figure 7c ) . 



* is the Kleene star denoting finite repetitions of a regular expression (see for 
example [HMU03J for details). Finally, for a configuration b we denote 6[ ni0 o[ the 
infinite subsequence of b starting from index ntocxi. The study presented in this 
paper leads to the result: 

Theorem 3. There exists an n in O(logiV) such that 

n(N) [nM e (( + ^ P ----2-l\ (T 

Some additional work will prove the more precise statement (named Snowball 
Conjecture in [PRllbJ): 

Theorem 4. There exists an n in O(logiV) such that 

7r(iV) [ri)Oo[ e(p-...-2-l)*[0](p-...-2-l)*0 w 
where [0] stands for at most one symbol 0. 



Graphical representations of the Theorems are given on figures 7b and 7c 



The example of 7r(2000) for p = 4 is given on figure 9c of appendix O 



S 



2 Once upon a time, a snowball 



We consider the parameter p to be fixed. The proof of Theorem [3] works as follows. 
We begin by establishing a relation between the height difference representation 
and the shot vector, leading to the construction of a DDS in Z p+1 , such that the 
orbit of a well chosen point according to the number of grains N describes the 
fixed point ir(N) (subsection |2.1[ ). The aim is then to prove the convergence of 
this orbit in O(logiV) steps, such that the shot vector values it embeds involve 
the form described in Theorem [3] (subsections 2.2 2.3 and 2.4). 



2.1 Linking height differences and shot vector in a dynamical 
system 

In this subsection we construct a discrete dynamical system in Z p+1 such that 
the orbit of a particular point (chosen according to N) describes tt(N). 

A fixed point n(N) can be represented as a sequence of height differences 
(bi)ie\n (ie., bi = Tr(N)i for all i) or as a shot vector (ai)i e ^. Those two repre- 
sentations are obviously linked in various ways. In particular for any n we can 
compute the height difference at index n provided the number of firings at n—p, 
n and n+1 because b n is initially equal to (the case n = is discussed below) , 
and: a firing at Ti — p increases h n by 1; a, firing at ti decreases b n by p H- lj a 
firing at n + 1 increases bn by p. Therefore, bn 

= a n - P — (p + l)a n +pa n+ \, with 
< b n < p since ir(N) is a fixed point, and thus 

p p p 

This equation expresses the value of the shot vector at position n + 1 according 
to its values at positions n—p and n, and a bounded perturbation < ^ ^ 1. 
We take as an initial condition a_ p = N and cij = for — p < i < 0, representing 
the fact that the column is the only one receiving N units of height difference. 

Remark 5. Note that a n+ \ e N, thus —a n - p + (p + l)a n + b n = mod p. As 
a consequence, the value of b n is nearly determined: given a n - p and a n , there 
is only one possible value of b n , except when — a n _ p + (p + l)a n = mod p in 
which case b n equals or p. 

For example, consider 7r(2000) for p = 4 (see appendix [C]) . We have a§ = 120 
and ci4 = 189, so —04 + 5as = 411 = 3 mod p. From this knowledge, 6g is 
determined to be equal to 1, so that ag = — ^04 + fag + = 103 is an integer. 

We rewrite this relation as a linear system we can manipulate easily. a n+ \ is 
expressed in terms of a n _ p and a n , so we construct a sequence of vectors (A^)^ 
with X{ 6 N p+1 and such that X n = (a n _ p , a n _ p+ i, . . . , a n ) T where u 7 " stands 
for the transpose of v. Note that we consider only finite configurations, so there 
always exists an integer no such that X n = for no ^ n, with = (0, . . . , 0) T . 
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Given X n and b n we can compute X n+ \ with the relation 





( ° 


1 





o \ 




(o\ 


X n+ i = AX n + — J with A = 
P 










1 






1 


J = 








\-i 








p+1 




w 



in the canonical base B = (ep, ei, . . . , e p ), with A a p + 1 square matrix^] 

This system expresses the shot vector around position n + 1 (via X n+ \) in 
terms of the shot vector around position n (via X n ) and the height difference at 
n (via b n ). Thus the orbit of the point Xq = (N, 0, . . . , 0, do) in N p+1 describes 
the fixed point tt(N). 

Note that it may look odd to study the sequence (6j)j £ iM using a DDS that 
presupposes a knowledge on (&i)j<=[M- It is actually helpful because of the under- 
lined fact that the value b n is nearly determined (remark [5]): in the following we 
will make no supposition on the sequence (&i)j 6 tM (except that ^ hi ^ p for 
all i) and prove that the system converges exponentially quickly in N, such that 
from an n in C(log N) the sequence (6j)j^ n is determined to have a regular wavy 
shape. 

The system we get is a linear map plus a perturbation induced by the discrete- 
ness of the values of height differences. Though the perturbation is bounded by a 
global constant at each step (b n ^ p for all n since tt(N) is a fixed point), it seems 
that the non-linearity prevents classical methods to be conclusive concerning the 
convergence of this model. 

We denote (j> the corresponding transformation from Z p+1 to Z p+1 , which 
is composed of two parts: a matrix and a perturbation. Let R(x) = + 



-x* 



, <r , , the characteristic polynomial of A is (1 — or) R(x). We 
can first notice that 1 is a double eigenvalue. A second remark, which helps to 
get a clear picture of the system, is that all the other eigenvalues are distinct and 
lesser than 1 from Lemma [9] (using a bound by Enerstrom and Kakeya |GG97| ). 
Therefore there exists a basis such that the matrix of <j) is in Jordan normal 
form with a Jordan block of size 2. Then, we could project on the p — 1 other 
components to get a diagonal matrix for the transformation, hopefully exhibiting 
an understandably contracting behavior. We won't exactly follow this approach, 
but will use these remarks in subsection 12.31 

We tried to express the transformation (j) i n a basis such that its matrix 
is in Jordan normal form, but we did not manage to handle the effect of the 
perturbation expressed in such a basis. Therefore, we rather express in a basis 
such that the matrix and the perturbation act harmoniously. The proof of the 
main Theorem is done in three steps: 

1. the construction of a new dynamical system: we first express is a new 



basis B' , and then project along one component (subsection 2.2); 



3 As a convention, blank spaces are 0s and dotted spaces are filled by the sequence induced 
by its endpoints. 
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2. the behavior of this new dynamical system is easily tractable, and we will 
see that it converges exponentially quickly (in 0(log N) steps) to a constant 



vector (subsection 2.3); 



3. finally, we prove that as soon as the vector is constant, then the wavy shape 
of Theorem [3] takes place (subsection 2.4). 



2.2 A new version of the dynamical system . . . 

From the dynamical system X n +\ = AX n + in the canonical basis B, we 
construct a new dynamical system for (ft in two steps: first we change the basis 
of Z p+1 in which we express (ft, from the canonical one B, to a well chosen 
B'; then we project the transformation along the first component of B' . The 
resulting system on Z p , called averaging system, is very easily understandable, 
very intuitive, and the proof of its convergence to a constant vector can then be 
concluded straightforwardly. Let 



B' 



(I 
1 1 

V 1 



0\ 




B'~ 



(I 



o\ 



-1 I) 



be square matricea^Jof size p + 1. B' = (eg, . . . , e' p ) (with e[ the (i + l) th column 
of the matrix B') is a basis of Z p+1 , and we have 



B'- 1 X,, 



n+l 



B'- 1 AB' B'- l X n 



Kb 1 - 1 j 

v 



with 



X' n — B' 1 X n 



A' = B' AB' 



A 1 







J' 



B'- l J 



M 


VV 



We now proceed to the second step by projection along e' . Let P denote 
the projection in Z p+1 along e' onto {0} x Z p . We can notice that e' is an 
eigenvector of M, hence projecting along e' Q simply corresponds to erasing the 
first coordinate of X[. For convenience, we do not write the zero component of 
objects in {0} x Z p . 

The new DDS we now have to study, which we call averaging system, is 



Y n+l =MY n 



-K 



P 



(1) 



4 As a convention, blank spaces are 0s and dotted spaces are filled by the sequence induced 
by its extremities. 
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with the following elements in Z p (in {0} x IP) 

/0 1 0\ 



Y n — P X' 



M = PA' 



1 

V 1 1 ••• -) 

\p V V' 



K = PJ' 





W 



Let us look in more details at Y n and what it represents concerning the shot 
vector. We have X„ - ■ ■■ ■■■ -■ ,,! 



&n—pi 0"n— p+lj • • • j On J ; thus 

^Q-n— p+1 Q"n—p^ 



Y n = PB'- 1 X n 



O-n-l — O-n-2 

\ a n — a n -\ ) 



and for initialization Yn 



(-N\ 




V a J 



Y n represents differences of shot vector, which may of course be negative. In the 
next subsection |2.3| we will see that the averaging system is easily tractable and 
converges exponentially to a constant vector. Subsection |2.4| then concentrates 
on implications following the presence of a constant vector i.e., the creation and 
maintenance of a wavy shape. 



2.3 



which is easy to tackle 



The averaging system is understandable in simple terms. From a vector Y n of 
IP, Y n+ i is obtained by 

1. shifting all the values one row upward; 

2. for the bottom component, computing the mean of values of Y n , and adding 
a small perturbation (a multiple of jjj between and 1) to it. 

Remark 6. (li)ieiN are once more integer vectors, hence the perturbation added 
to the last component is again nearly determined: let y n denote the mean of the 
values of Y n , we have (y n |y)eZ and < ^ ^ 1. Consequently, if y n is not 
an integer then b n is determined and equals p(|y n l — Vn), otherwise b n equals 
or p. 

For example, consider 7r(2000) for p = 4 (see figure 9a of appendix |Cj be care- 
ful that it pictures a n — a n+ \ at position n). We have Y\s = (—3, —5, —7, — 7) T , 
then ?/i3 = — ^ and 613 is forced to be equal to 2 so that Y14 = (—5, —7, —7, — 5) T 
is an integer vector. 

We can foresee what happens as we iterate this dynamical system and new 
values are computed: on a large scale — when values are large compared to p — 
the system evolves roughly toward the mean of values of the initial vector Yn, 
and on a small scale — when values are small compared to p — the perturbation 
lets the vector wander a little around. Previous developments where intending 
to allow a simple argument to prove that this tiny wanderings does not prevent 
the exponential convergence towards a constant vector. 

The study of the convergence of the averaging system works in three steps: 



12 



(i) . state a linear convergence of the whole system; then express Y n in terms of 

Yq and (bi)o^ n ; 

(ii) . isolate the perturbations induced by (&i)o^i^n and bound it by a constant; 

(hi) . prove that the other part (corresponding to the linear map M) is contract- 
ing exponentially quickly. 

From (ii) and (iii), a point evolves exponentially quickly to a ball of constant 
radius, then from (i) this point needs a constant number of more iterations in 
order to reach the center of the ball, that is, a constant vector. 

For the sake of clarity, proofs of Lemmas are only sketched within the proof 
of the following Proposition, and are fully proven just below. 

Proposition 7. There exists an n in 0(\ogN) such that Y n is a constant vector. 

Proof. Let y n (respectively y n , y n ) denote the mean (respectively maximal, min- 
imal) of values of Y n . We will prove that y n — y^ converges exponentially quickly 
to 0, which proves the result. 

We start with Y Q = (-N, 0, . . . , 0, a ) T , thus y - y Q = N + a 2±I TV since 

ao ^ — (recall that ao is the number of times column has been fired). 

This proof is composed of two parts. Firstly, the system converges exponen- 
tially quickly on a large scale. Indeed, when y n — y^ is large compared to p, it is 
easy to be convinced that shifting the values and padding with the mean value 
plus a small perturbation decreases exponentially the gap between the values 
(because the perturbation is negligible). In other words, the value of y n — y n 
decreases exponentially to a ball of constant radius: 



Lemma 8. There exists a constant a and a no in C?(log N) s.t. y^—y < a. 

Proof sketch. Complete proof below. Since Y n converges towards the mean 
of its values, we consider the evolution of the distance to the mean vector 
associated to Y n . We express Y n in terms of Yq and (&i)o=si^m anci study 
separately the parts involving respectively Yq (part 1), and (6j)o^i^n (part 2). 
We first prove that part 1 converges exponentially quickly to since it consists 
in multiplying Yq by the n th power of a matrix which eigenvalues are the set 
of roots of R(x) (see Lemma [9| and 0, thus a contracting matrix. Then we 
easily upper bound part 2 by a constant a independent of N, which completes 
the proof. □ 



Lemma 9. Let R(x) = x^ 1 + ^x^ 2 + • • • + |x + |. 
R{x) has p — 1 distinct roots Ai, . . . , A p _i and for all i, Aj ^ 

Proof sketch. Complete proof below. It uses a bound by Enestrom and Kakeya 
(see for example [GG97j), and the fact that R{x) and R'{x) are coprime. □ 

Secondly, on a small scale, we use the fact that the system converges linearly 
(values of Y n are integers). 
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Lemma 10. The value of y n — y n decreases linearly: if y y n , then there 
is an integer c, with < c < p such that y n+c — U n+C < V n ~ U n - 

Proof sketch. Complete proof below. The result follows the observation of the 
strict inequalities y^ < y n < y n . y n is used for the computation of the bottom 
component of Y n+ i, hence the idea is that all the values of Y n+p (recall the 
shifting of values) lies strictly in between the values of Y n . □ 

To conclude, we start with y — y Q in 0(N), we have a constant a and a uq in 
O(logiV) such that y nQ — y^ < a thanks to the exponential decrease on a large 
scale (Lemma |8|. Then after p iterations the value of y n+p — D n+p is decreased 



by at least 1 (Lemma 10), hence there exists (3 with (3 < pa such that after j3 
more iterations we have y no+ p — y n + ^ = 0. Thus Y no+ /3 is a constant vector, 
and n + (3 is in C(logiV). □ 

In this proof, neither the discrete nor the continuous studies are conclusive by 
themselves. One one hand, the discrete study gives a linear convergence but not 
an exponential convergence. On the other hand, the continuous study gives an 
exponentially convergence towards a constant vector, but in itself the continuous 
part never reaches the constant vector but tends asymptotically towards it. It is 
the simultaneous study of those modalities (discrete and continuous) that allows 
to conclude. 

Remark 11. Note that for p = 1, the averaging system has a trivial dynamic. For 
p = 2, the behavior is a bit more complex, but major simplifications are found: 
the computed value is equal to the mean of 2 values, hence the difference y n — y n 
decreases by a factor of 2 at each time step. The key simplification arising for 
p = 2 is the fact that this decrease arises straightforwardly at each time step. 

Lemma 8. There exists a constant a and a no in 0(logN) s.t. y nQ — < a. 

Proof. We start with Y = (-N, 0, . . . , 0, a ) T , thus y - y Q = N + a «S 
The relation linking Y n to Y n+ \ is 

Y n+l = MY n + b -^K 
V 

Since we want to prove that Y n converges to a constant vector close to the 
mean of its values, we want to consider the evolution of the distance to the mean 
vector associated to Y n via the following application D: 



with y n = - V y n 
p ^r 1 

The aim is thus to prove that there exists an no in O(logiV) such that the 
norm of DY no is bounded by a constant. 









( y no -y n \ 


Y n = 










\Vnp-J 




\2/n p _i ~ Vn) 
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We can notice that Y n — DY n = (y n , . . . ,y n ) T > therefore M (Y n — DY n ) = 
(Y n -DY n ), and we have (D MY n ) - (D M DY n ) = 0. Consequently, we can 
consider the relation 

D Y n+1 = DM DY n + — DK 

V 

and express Y n in terms of Yq and (6j)o^i^n : 

-. 71—1 

DY n = (DM) n DY + - 2 HDMY'-^DK 

p i=0 

In order to prove the result, we will see that the linear map Z *—> (D M) Z is 
eventually contracting, hence it converges exponentially quickly to 0, its unique 
fixed point ( |KH96| Corollary 2.6.13). That is, (DM) n DY converges to 
exponentially quickly. It then remains to upper bound the norm of the remaining 
sum by a to get the result. 

To prove that the map Z *—> (D M) Z is eventually contracting, it is enough 
to prove that its spectral radiu^ is smaller than 1 ( |KH96j Corollary 3.3.5). 
This part is detailed in Lemma [12] below, using the fact that M is a companion 
matrix which eigenvalues are simply upper bounded with a result by Enestrom 
and Kakeya |GG97| . 

Since y — y Q is in 0(N), ||Z) Yolloo is also in 0{N) and there exists an uq in 
O(logiV) suchlhat || (D M) n ° D Yb||ao < 1. 

It remains to upper bound the summation by a constant (we recall that for 
a matrix A, \\A\lca = sup \\A x||co for ||x||co = 1): 

i n £ biiDM^-^DK 

i=0 



for some constant f3 independent of N. Finally, we have 

-. n— 1 

\\DY no \\ x < \\(DM) no DY \\ x + \\-Y b i {DM) nr ~ 1 ~~ i DKU < f3 

and the fact that y no ~ U n) ^ 2||Z)y no || co completes the proof with a = 2/3. □ 

Lemma 9. Let R(x) = x p ~ l + ^x p ~ 2 H ^\ x +\- 

R(x) has p — 1 distinct roots Ai, . . . , A p _i and for all i, Aj ^ 

Proof. The distinctness of the roots of S(x) = px p ~ l R(^) = + 2x p ~ 2 + • • • + 
{p — l)x + p implies the distinctness of the roots of R(x). The distinctness of 
the roots of S(x) comes from the fact that S(x) and S'(x) are co-prime. With 
a = ^Tj x + 3 and b = ^+T) x2 ~ ^+T) X ^ we g et aS{x) + bS'{x) = 1. Therefore 
from Bezout GCD(R{x),R'{x)) = 1, which implies the result. 

the maximal absolute value of an eigenvalue of DM. 



oo 



n -l 

sc i £ p\\DM\ 



no— 1— i 



i=0 
1 



oo 



DK\ 



00 



1-pMlloo 

B-l 



\\DK\ 



00 
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For the second part of the lemma, a classical result due to Enestrom and 
Kakeya (see for example [GG97J) concerning the bounds of the moduli of the 
zeros of polynomials having positive real coefficients states that all the complex 
roots of R(x) have a moduli smaller or equal to 2^1. □ 

Lemma 10. The value ofy n — y n decreases linearly: if y # y n , then there is 
an integer c, with < c < p such that y n+c — U n+C < V n ~ V n - 

Proof. Let Y n = (Y no , . . . , Y rip _ 1 ). If y n y n , that is, if the vector is not constant, 
the mean value is strictly between the greatest and smallest values: y n <y n < y n . 
Consequently y n < {Y n+ i p _ 1 = y n + ^ < y n . This reasoning applies while 
y n+i 7^ y n +i- If we consider p iterations, we either have y n+c = U n+C f° r some 
c ^ p; or, we use the facts that y < V n+ i ano - Vn ^ Vn+l> therefore we have 
y < Y n+Pi sC y n for all i (because shifting the values leads to Y n+Pi = Y n+ \ + i p _ 1 
for all «S i < p - 1, and y n V n+1+i < Y n¥l+ip _ 1 «S Vn+i+i < V n ), thus 

y n < m„ +p and y« +P < i/n- □ 

Lemma 12. T/ie spectral radius of DM is strictly smaller than 1. 

Proof. M is a compagnon matrix, its characteristic polynomial is 

p-l 



2 ~x fc = (x- l)R(x) 



k= Q P 

with i?(x) = + 2--^x p ~ 2 + • — h + i. From Lemma j9j we know that R(x) 

has p — 1 distinct roots Ai, . . . , A p _i, all comprised between ^ and The set 
of eigenvalues of M is thus M\ = {1, Ai, . . . , A p _i}. We will prove that the set of 
eigenvalues of the matrix DM is DM\ = {0, Ai, . . . , A p _i}. 

Let vo, . . . , Vp-% be non null eigenvectors respectively associated to the eigen- 
values 1, Ai, . . . , A p _i of M. The case v\ is particular and allows to conclude that 
is an eigenvalue of DM. The other eigenvectors of M lead to the conclusion 
that DM also admits the eigenvalues Ai, . . . , \ p -\. 

• DM vq = Dvq since the associated eigenvalue is 1, and Dvq = because 
the eigenspace associated to the eigenvalue 1 is the hyper plan of constant 
vectors. As a consequence, is an eigenvalue of DM. 

• For the other eigenvectors, that is, for 1 < i < p — 1, let C, be the con- 
stant vector with all its components equal to - Xlfc=o Vi k > with v i k the k th 
component of the vector v j. Vi — c\ # 0, and 

DM [vi -Ci) = D (M v,i — M a) 
= D (Ai vi - a) 
= Ai D Vi — D Ci 
= Xi(vi - Ci) - 

where the last equality is obtained from the fact that by definition of D we 
have D Vi = Vi — Ci. As a consequence, Ai is an eigenvalue of DM. 

Finally, DM\ = {0, Ai, . . . , X p — 1} and the spectral radius of DM is smaller 
or equal to □ 
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2.4 . . . and allows to conclude 

Lemma [7] shows that there exists an n in O(logiV) such that Y n is a constant 
vector. In this subsection, we prove that if Y n is a constant vector, then the 
shape of the sand pile configuration is wavy from the index n, exactly as it is 



stated in the main Theorem and on figure 7b 



Lemma 13. Y n is a constant vector of Z p implies 

Proof. From this index n for which Y n is a constant vector, we will see that the 
sequence (6j)j^ n is determined, or more accurately nearly determined. The proof 
works as follows. We will first see that if Y n is a constant vector, then the value 
of b n is or p. if it is then Y n+ i is again a constant vector; if it is p, then 
the sequence (bi) n <zi <n+p is determined to be equal to (p,p — 1, . . . , 1). Lemma 
[j~6l states that there is no sequence of more than p + 1 values in the sequence 
(6j)j e jsl, hence we may have a sequence of at most p + 1 values (during which 
Yi remains constant), eventually followed by a value p which triggers the pattern 
(p,p— 1, . . . , 1). Finally, at the end of the sequence (p,p— 1, . . . , 1), we are back 
to a constant vector Yi, and the reasoning can be repeated until the end of the 
configuration (the W ), hence the result. 

We concentrate on the sequence of values of Y{. The fact that its components 
are integers, and especially the last one, will play a crucial role in the determi- 
nation of the value of hi because ^ b% p (let us recall that the sequence 
b{ is the height difference representation of the fixed point with N grains, i.e., 

We start from the hypothesis that Y n = (a, . . . , a) T , thus from the averaging 



b 

IX T 

vector and a is an mleger, hence b n equals or p. 



system's equation (1) we have Y n+ \ = (a, . . .,a,a + -y) . Y n+ \ is an integer 



If b n = then Y n+ \ = (a, . . . , a) T and we are back to the same situation, 
the dilemma goes on: the value of 6 n +i is not determined, it can be or p. 

If b n = p then Y n+k+ i = (a, . . . , a, a + 1) T fr om the relation above. A 
regular pattern then emerges: 

- ifY n+ i = (a,...,a,a + l) T ,thenY n+2 = (a, . . . , a, a + 1, pa+1 + b " +1 ) T 
and it determines b n+ \ = p — 1 so that Y n+ 2 = (a, . . . , a, a. + 1, a + 1) 
is an integer vector; 

- if Y n+2 = (a, . . . ,a, a+1, a+l) T , then Y n+2 = (a, . . . , a,a+l, P a+2 + ! 
and it determines b n+ 2 = p — 2 so that Y n+ z = (a, . . . , a, a + 1, a + 
1, a + 1) T is an integer vector; 

- et cetera we have b n+ i = p — i for ^ i < p, and eventually Y n+p = 
(a + 1, . . . , a + 1) T is a constant vector (note that Yq has a negative 
mean, hence a is negative, which is consistent with the a + 1 we 
obtain). 



n + 2 ' 
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Let us illustrate the developments above on the following picture, where ar- 
rows are labeled by values of b%. 



p-X 



o— o- 

i 



When Y n is constant we are in the grey node, then it is either possible that b n = 
in which case we are back in a situation where Y n +\ is constant, or b n = p in 
which case b n+ \ = p — 1, . . . , b n+p -i = 1 and we are back in a situation where 
Y n+p is a constant vector. The last argument is that from Lemma 16 it is not 
possible to have a sequence of more than p + 1 consecutive 0, hence it is not 
possible to loop more than p + 1 times on the grey node. 

To conclude, if Y n is a constant vector we are on the grey node, and then the 
set of possibilities for the sequence (6j) n ^j is described by paths on the automata 
above. Furthermore it is not possible to loop more than p + 1 consecutive times 
on the grey node, consequently the set of possibilities for the sequence (bi) n <zi 
is exactly the one described in the statement of the Theorem. Note that it is 
eventually possible to loop infinitely on the label 0, when there is no more grains 



on columns since Lemma 16 won't apply anymore, hence the W . □ 



2.5 Proofs of Theorem [3] and Theorem S] 

It remains to put the two pieces together to prove Theorem [3] 
Theorem 3. There exists an n in 0(logN) such that 




7r(N) [nM e -\-0 l -p-...-2-l) (f 



Proof. From Proposition [7J there exists an n in 0(logiV) such that Y n is a 
constant vector. We can therefore apply Lemma [13] for this n and the result 
follows. □ 

An example of the various representations of fixed points (X n , Y n and wavy 
shape) for 7r(2000) with p = 4 is given in appendix [U] 

In order to prove the refinement of Theorem [3] it seems necessary to overcome 
the "static" study — direct, for any fixed point — presented here and consider 
the dynamic of sand grains on columns beyond 0(}ogN), from 7r(0) to tt(N). 
We recall that ir(N) can be computed inductively, using the relation 

for all k > 0, 7r(vr(A; - = ir(k) 

where denotes the configuration obtained by adding one grain on column 
of a. We start from 7r(0) and inductively compute 7r(l),7r(2), . . . , tt(N— 1) and 
tt(N) by repeating the addition of one grain on column 0. The sequence of firings 



from 7r(fc — 1)^° to ir(k) is called the k th avalanche (see subsection 1.4). 

An interesting Corollary of Theorem [3j proven below, is that for any p and 
N, there is an index n in C(log N) such that the N th avalanche fires a set of con- 
secutive columns (without missing any) on the right of that index n (conjectured 
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in |PRlla] and proven only for p = 2). We now formally introduce this impor- 
tant property on avalanches, which, when it is verified starting from an index 
n, leads to regularities in the avalanche process beyond column n (see |PRlla] 
and |PRllb] for details). We will see that a Corollary of Theorem [3] is that this 
property is verified on the k th avalanche starting from an index in 0(logA;). 

We say that there is a hole at position i in an avalanche s k if and only if i ^ s k 
and (i + 1) £ s k . An interesting property of an avalanche is the absence of hole 
from an index /, which tells that, 

there exists an m such that [ for all i with I < i ^ m, we have i e s k 

for all i with m < i, we have i £ s k 

namely, from column I, a set of consecutive columns is fired, and nothing else. 
We say that an avalanche s k is dense starting from an index I when s k contains 
no hole i with i ^ I. We have already explained in |PRllaj that this property 
induces a kind of "pseudo linearity" on avalanches, that it somehow "breaks" the 
criticality of avalanche's behavior and let them flow smoothly along the sand pile. 
Let us introduce a formal definition: 

Definition 14. C'(p, k) is the minimal column such that the k th avalanche is 
dense starting at £' (p, k) : 

£'(p, k) = min{/ e N | 3m e N such that V/ ^ i ^ m, i e s k and Mi > m, i ^ s k } 

Then, the global density column C(p, N) is defined as: 

C(p,N) = max{£ (p,k) | k ^ N} 

The global density column C(p, N) is the smallest column number starting 
from which the N first avalanches are dense (contain no hole). A Corollary of 
Theorem [3] is: 

Corollary 15. For all parameter p, C(p,N) is in (D(logN). 

Proof. Let p be a fixed parameter. In order to prove the result, we will prove 
that if 



vr(% nM e +0* -p-...-2-l 




then the next avalanche s k+l is dense starting from n + c with c i-£-, 
C'(p, k + 1) ^ n + c. Since such an n is in 0(log k) for all k from Theorem |3j the 
Corollary follows. 

Let us consider the following case disjunction. 

• If the k+l th avalanche ends before column n, formally if maxs fe+1 < n, 
then obviously C'(p, k + 1) < n because for all /, C'(p, I) ^ maxs'. 

• If the k + l th avalanche ends beyond column n, formally if maxs fc+1 ^ n, 
then we will consider the dynamic of the avalanche process s k+l on columns 
greater than n. There is an intuitive notion of time during an avalanche, 
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induced by the order of the firings: s k+1 is the fired column at time t. For 
the sake of clarity, let us precise our hypothesis given by Theorem [3] on the 
shape of 7r(fc)r Tl)00 [. Let I and z\, . . . , zi be the integers such that 

21 Z 2 Zl 

A k )[nM = 0~ ... • 0-p-. . .-2-1-0 0-p-. . .-2-1 • . . . • 0-p-. . .-2-l-0 w 

The aim is to prove that s k+1 is dense starting from n + z±, which is proven 
by induction on the following Statement: 

Statement. For all i 5= n + z\, if i $ s k+1 then (i + l) $ s k+1 . 

Let us prove this Statement for alH > n + Z\, by considering two cases. 

— If 7r(k)i > 0, we prove the contraposition of the Statement: if (i + l) e 
g k+i n ien i e s k+l . The firing at (i+l) gives p units of height difference 
to column i, which becomes unstable and thus i e s k+1 . 

— If n(k)i = 0, we prove the Statement by contradiction: let us suppose 
that (i + l) £ s k+1 . Column i can: receive p units of height difference 
from (i+l), and 1 unit of height difference from (i—p). By hypothesis, 
i receives p units of height difference thus (i—p) ^ s k+l otherwise i 
receives p+1 units of height difference during the avalanche process 
and i e s k+1 . 

According to the shape of 7r(k) beyond column n, there are two possi- 
bilities for the value of Tr(k)i+\. Let us show that in any case we have 
(i-p+1) e s k+1 . 

* If 7r(/c)j+i = 0, then since (i+l) is fired it should receive at least p+1 
units of height difference, which is possible only if (i— p+1) £ s k+1 . 

* If 7r(A;)j + i = p, then we prove it in two steps. Firstly, from the 
locality of the rule (it involves columns at distance at most p) and 
the stability of ir(k), a column can't be fired if none of the p pre- 
ceding columns has been fired. Secondly, for all j with < j < p, 
we have vr(fc)j + i_j < p thus by induction on fact that such a 
column can't be fired if its successor (the consecutive column on 
its right) has been fired, i + l is fired before any of the p preced- 
ing columns. Consequently, it is not possible to fire any column 
greater than i + l before column i + l itself. We can therefore con- 
clude that i + l received units of height difference from the left of 
the configuration, that is, from (i— p+1). 

We proved that (i— p + 1) e s k+1 . Now, if (i + l) e s k+1 then the 
z L associated to column i is at most p—1, otherwise the avalanche 
process can't go beyond this flat surface (again, from the locality of 
the rule). As a consequence, 7r(/c)j_ p > 0, and since (i— p + 1) £ s k+1 
it gives p units of height difference to i—p thus it becomes unstable, 
but (i—p) ^ s k+1 1 which is a contradiction. 

□ 
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Theorem 4. There exists an n in 0(logN) such that 



7r(iV) [niOo[ e(p-....2-l)*[0](p-...-2-l)*0 w 
where [0] stands for at most one symbol 0. 

Proof. We prove a slightly stronger result, that there is an n in O(logiV) such 
that 

n(N) [nM e M- 0M (p • . . . • 2 • l)*[0](p • . . . • 2 • 1)*0" 

The proof works by induction on N, with the base case coming from Theorem 
[31 For the induction, suppose the result holds for ir(N) with the wavy pattern 
l (p- . . . ■ l) x 0(p- . . . • l) v , with i ^ 0, x > and y > 0, beginning at n in 0(log N). 
In addition, there is a minimal index 2 such that the (N+l) th avalanche fires a set 



of consecutive columns, from I to m, and from Corollary 15 this / is in O(logiV). 
Without loss of generality, we consider that I < n — p. 

• if m < n + i, either the wavy pattern is not affected (if m =S n — p), or 
Theorem |3] ensures that there will always exist an index n in C(log N) such 
that the wavy pattern is preserved; 

• if n + i < m, we first note that i =fi p + 1 because an avalanche cannot 
"cross" a set of p + 1 consecutive 0. In this case, the avalanche triggers a 
chain reaction that stops on the last wave before the is encountered, and 
thus spans the x first waves, i.e., m = n + {x — l)p. A simple look at the 
transition rule indicates that the firing of columns n — p to m shifts the 
wavy pattern to the following: 0*(p • . . . • l) x ^ 1 0(p ■ . . . ■ l) y+1 . 

In any case, the result holds for tt(N + 1). □ 



3 Conclusive discussion 

The result of this paper stresses the fact that sand pile models are at the edge 
between discrete and continuous systems. Indeed, when there are very few sand 
grains, each one seems to contribute greatly to the global shape of the config- 
uration. However, when the number of grains is very large, a particular sand 
grain seems to have no importance to describe the shape of a configuration. The 
result can be interpreted has a separation of the discrete and continuous parts 
of the system, despite it is not clearly stated. On one hand, the left seemingly 
unordered part, interpreted as reflecting the discrete behavior, prevents regular- 
ities to emerge. On the other hand, the right and ordered part, interpreted as 
reflecting the continuous behavior, lets a regular and smooth pattern to come 
into view. 

Nevertheless, the separation between discrete and continuous behaviors may 
be challenged because the continuous part emerges from the discrete part. We 
have two remarks about this latter fact. Firstly, the consequence seems to be a 
slight bias appearing on the continuous part: it is not fully homogeneous — that 
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is, with the exact same height difference at each index — which would have been 
expected for a continuous system, but a — very small — pattern is repeated. This 
remark actually seems to be misleading because this exact same bias disappears 
when we consider an initial configuration with p consecutive indices with N units 
of height difference, thus it looks like the bias comes from the gap between the 
unicity of the initial column compared to the parameter p of the rule. Secondly, if 
we consider the asymptotic form of a fixed point, the relative size of the discrete 
part is null. This, regarding the intuition described above that when the number 
of grains is very large then a particular grain has no importance, is satisfying. 

Let us save the last words to a distracting application to a famous paradox. 
Someone who has a very little amount of money is called poor. Someone poor 
who receives one cent remains poor. Nonetheless, if the increase by 1 cent is 
repeated a great number of time then the person becomes rich. The question is: 
when exactly does the person becomes rich? An answer may be that richness 
appears when money creates waves... 
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A There is no plateau of length larger than p + 1 



A plateau is a set of at least two non empty and consecutive columns of equal 
height. The length of a plateau is the number of columns composing it. 

Lemma 16. For any N and any configuration a such that (N, W ) — > a, in a 
there is no plateau of length strictly greater than p + 1. 

Proof. This proof proceeds by contradiction, using the fact that configurations 
are sequences of non-negative integers (Hi). Suppose there exists a configuration 
a reachable from (N, 0^) for some N, such that there is a plateau of length at 
least p + 2 in a. Since there is no plateau in the initial configuration (JV, W ), 
and there is a finite number of steps to reach a, there exists two configurations 
p and t such that p —* r and such that there is a plateau of length at least p + 2 
in r, and none in p (T^)- Let k be the leftmost column of the plateau of length 
at least p + 2 in r, i.e. for all column j between k and k + p, tj = (H3). We 

will now see that there is no i such that p —> r, which completes the proof. 

• ii i < k — poii>k+p + l then a firing at i has no influence on columns 
between k and k + p + 1 and there is a plateau of length at least p + 2 in 
p, contradicting %2- 

• iffe — p ^ i ^ k then according to the rule definition we have T{ +p = 
Pi+p-1 — 1) an d from %3 /9j +p = therefore Tj +P < which is not possible 
(recalled in Hi). 

• \lk<i^k+p+l then according to the rule definition we have Tj_i = 
Pi_i — p and from pi—\ = therefore Ti—\ < which again is not possible 



We give bounds for the maximal index of a non-empty column in the fixed point 
tt(N) according to the number iV of grains, denoted w(N). w(N) can be inter- 
preted as the support or width or size of tt(N). We consider a general model 
KSPM(p) with p a constant integer greater or equal to 1. A formal definition of 



Lemma 17. The support of 'ir(N) is in Q(y/N). 

Proof. The support of ir(N) is denoted w(N). 

Lower bound: ir(N) is a fixed point, therefore by definition for all index i we 
have 7r(A r )j ^ p. Then, 



from Hi. 



□ 



B 



The support of ir(N) is in 6(V^V) 




w(N) 

n < y] P ■ * = P 

i=0 



w{N) • (w(N) + 1) 
2 



< p 2 {w{N) + l) 2 
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Figure 8: The support of tt(N) is in 0(yjV). It is lower bounded by the fact that 
each height difference on a stable configuration is at most p, and upper bounded 
by the fact that there is no plateau of length larger than p + 1. 



hence \\fN - 1 < w(N). 

Upper bound: From Lemma 16 there is no plateau of length greater than 
p + 1. Therefore, for w(N) > p we have 



I I 

L p+1 J 



N > J] (P + 1) • * > (P + 1) 



p+1 



1 p+l 



> 



i=0 



w(JV) 

p + 1 



hence (p + 1)VN + p + 1 > io(jV). 



□ 



C tt(2000) for p = 4 

Figure [9] presents some representations of 7r(2000) for p = 4 used in the develop- 
ments of the paper: differences of shot vector (figure 9a), shot vector (figure 9b) 



and height (figure 9c). 



vr(2000) 



(4,0,4,1,3,2,4,1,1,3,4,3,4,2,0,1,4,2,2,1, 

4, 3, 2, 1,0, 4, 3, 2, 1,4, 3, 2, 1,4, 3, 2, 1,4, 3, 2, 1,0") 



We can notice on figure 9a that the shot vector differences contract towards 
some "steps" of length p, which corresponds to the statement of Lemma [7] that 
the vector Y n becomes constant exponentially quickly (note that this graphic 
plots the opposite of the values of the components of Y n ). The shot vector 
representation on figure 9b" corresponds to the values of the components of X n , 



which we did not manage to tackle with classical methods. Figure 9c pictures the 
sand pile configuration on which the wavy shape appears starting from column 
20. 
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